Direct Characterization of Quantum Dynamics 
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The characterization of quantum dynamics is a fundamental and central task in quantum mechanics. This 
task is typically addressed by quantum process tomography (QPT). Here we present an alternative "direct char- 
acterization of quantum dynamics" (DCQD) algorithm. In contrast to all known QPT methods, this algorithm 
relies on error-detection techniques and does not require any quantum state tomography. We illustrate that, by 
construction, the DCQD algorithm can be applied to the task of obtaining partial information about quantum 
dynamics. Furthermore, we argue that the DCQD algorithm is experimentally implementable in a variety of 
prominent quantum information processing systems, and show how it can be realized in photonic systems with 
present day technology. 
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The characterization and identification of quantum systems 
are among the central modern challenges of quantum physics, 
and play an especially fundamental role in quantum informa- 
tion science fH, and coherent control J2l- A task of gen- 
eral and crucial importance is the characterization of the dy- 
namics of a quantum system that has an unknown interac- 
tion with its embedding environment. Knowledge of this dy- 
namics is indispensable, e.g., for verifying the performance 
of an information-processing device, and for the design of 
decoherence-mitigation methods. Such characterization of 
quantum dynamics is possible via the method of standard 
quantum process tomography (SQPT) SQPT consists 

of preparing an ensemble of identical quantum systems in a 
member of a set of quantum states, followed by a reconstruc- 
tion of the dynamical process via quantum state tomography. 
We identify three main issues associated with the required 
physical resources in SQPT. (i) The number of ensemble mea- 
surements grows exponentially with the number of degrees of 
freedom of the system, (ii) Often it is not possible to pre- 
pare the complete set of required quantum input states, (iii) 
Information concerning the dynamical process is acquired in- 
directly via quantum state tomography, which results in an 
inherent redundancy of physical resources associated with the 
estimation of some superfluous parameters. To address (ii), 
the method of ancilla-assisted process tomography (AAPT) 
was proposed [0] . However, the number of measurements is 
the same in SQPT and (separable) AAPT 0]. 

Here we develop an algorithm for direct characterization 
of quantum dynamics (DCQD), which does not require quan- 
tum state tomography. The primary system is initially en- 
tangled with an ancillary system, before being subjected to 
the unknown dynamics. Complete information about the dy- 
namics is then obtained by performing a certain set of error- 
detecting measurements. We demonstrate that for character- 
izing a non-trace preserving quantum dynamical map on n 
qubits the number of required experimental configurations is 
reduced from 2 4 ", for SQPT and separable AAPT, to 2 2 " in 
DCQD. E.g., for a single qubit, we show that one can fully 



characterize the quantum dynamics by preparing one of four 
possible two-qubit entangled states, and a Bell-state measure- 
ment (BSM) at the output. This is illustrated in Fig.Q]and Ta- 
ble U We also discuss the experimental feasibility of DCQD 
in a variety of quantum information processing (QIP) systems, 
and show how it can be realized with linear optics. 

In principle, the number of required experimental configu- 
rations in the AAPT scheme can be reduced by utilizing non- 
separable (global) measurements for the required state tomog- 
raphy, such as mutually unbiased basis (MUB) measurements 
J5[|, or a generalized measurement |@]. In general, these types 
of measurements require many-body interactions which are 
not available experimentally. Here, we demonstrate that the 
DCQD algorithm requires only 0(n) single- and two-body 
operations per experimental configuration. 

We demonstrate the inherent applicability of the DCQD 
algorithm to the task of partial characterization of quantum 
dynamics in terms of coarse-grained quantities. Specifically, 
we demonstrate that for a two-level system undergoing a se- 
quence of amplitude and phase damping processes, the relax- 
ation time T\ and dephasing time T2 can be simultaneously 
determined in one ensemble measurement. 

Quantum Dynamics. — The evolution of a quantum system 
(open or closed) can, under natural assumptions, be expressed 
in terms of a completely positive quantum dynamical map £ , 
which can be represented as 11J 
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Here p is the system initial state, and the {E m } are a set of (er- 
ror) operator basis elements in the Hilbert-Schmidt space of 
linear operators acting on the system, satisfying Ti(EjEj) = 
dSij. The {xmn} are the matrix elements of the superopera- 
tor x, which encodes all the information about the dynamics, 
relative to the basis set {E m } (U. For an n qubit system, 
the number of independent matrix elements in x K 2 4 " for a 
non-trace-preserving map, or 2 4 ' 1 — 2 2n for a trace-preserving 
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FIG. 1 : Schematic of characterization algorithm for single qubit case, 
consisting of Bell state preparations, applying the unknown quantum 
map, £ , and Bell-state measurement (BSM). 

map. The matrix x is positive Hermitian, and Trx < 1. Thus 
X can be thought of as a density matrix in Hilbert-Schmidt 
space, whence we often refer, below, to its diagonal and off- 
diagonal elements as "quantum dynamical population" and 
"coherence", respectively. 

Characterization of Quantum Dynamical Population. — 
Here we demonstrate how to determine {xmm} in a single 
(ensemble) measurement for a single qubit. Let us maxi- 
mally entangle two qubits A (primary) and B (ancillary) as 
|V') = (|0a0_b) + \IaIb))/\/^, and then subject only qubit 
A to a map £ . From here on, for simplicity, we denote £ (g> X 
by £. In this case the error basis {-E m }^ =0 becomes the iden- 
tity operator and the Pauli operators: {/, X, Y, Z}. The state 
\tp) is a +1 eigenstate of the commuting operators Z A Z B and 
X A X B , i.e., it is stabilized under the action of these "stabi- 
lizer operators", and it is referred to as a "stabilizer state" 
Any non-trivial operator acting on the state of the first qubit 
anticommutes with at least one of Z A Z B and X A X B , and 
therefore by measuring these operators we can detect an arbi- 
trary error on the first qubit, i.e., by finding the eigenvalues of 
either one or both operators to be —1. Measuring the observ- 
ables Z A Z B and X A X B is equivalent to a BSM, and can be 
represented by the four projection operators Pq — (0 + 1, 

Pi = P2 = |V>-)(V>~|. Pz = l<r>(0~|. where 

|</>±) = (|00) ± \11))/V2, |^±) = (|10) ± |01»/V2, form 
the Bell basis of the two-qubit system. The probabilities of 
obtaining the no error outcome /, bit flip error X A , phase 
flip error Z A , and both phase flip and bit flip errors Y A , 
on the first qubit, become p rn = Ti[P m £ (p)] = Xmm, for 
in = 0, 1, 2, 3 respectively. Therefore, we can determine the 
quantum dynamical population, {Xmm}m = Q, in a single en- 
semble measurement (e.g., by simultaneously measuring the 
operators Z A Z B and X A X B ) on multiple copies of the state 

\i>)- 

Characterization of Quantum Dynamical Coherence. — In 
order to preserve the coherence (Xm^n) in the quantum dy- 
namical process, we perform a set of measurements such that 
we always obtain partial information about the nature of the 
errors. The simplest case is a measurement of the projec- 
tion operators P+i = \(f> + )((f> + \ + \<p-)((j>-\ and P_i = 
\tp + )(ip + \ + \"4>~) | corresponding to the eigenvalues +1 
and —1 of measuring the stabilizer Z A Z B . The outcomes of 
this measurement represent the probabilities of no bit flip er- 
ror and bit flip error on qubit A, without telling us anything 
about a phase flip error. Therefore, we preserve only the co- 
herence between operators / and Z A , and also between the 
X A and Y A (which are represented by the off-diagonal el- 



ements X03 a nd X12, respectively). The required input state 
is a non-maximally entangled state \4>c) = a\00) + /3|11), 
with \a\ ^ |/?| 7^ 0, whose sole stabilizer is Z A Z B . Now, by 
separating real and imaginary parts, the probabilities of no bit 
flip error and bit flip error events become: Tr[P + i£ (p)] = 
Xoo + X33 + 2Re( X03 )(Z A ) and Ti[P^£(p)] = X n + 
X22 + 2lm(xi 2 ){Z A ), where {Z A ) = Ti(pZ A ) ^ (be- 
cause |a| 7^ |/3| 7^ 0), with p = \<fic) (4>c\- We already 
know the {xmm} from the population measurement described 
above, so we can determine Re(xo3) and Im(xi2). Af- 
ter measuring Z A Z B the system is in either of the states 
p±i = P±\£ (p)P±i/Tr[P±i£(p)]. Now we measure the 
expectation value of a (commuting) normalizer operator U 
(such as X A X B ) 0]. Thus we obtain Tr[f7p + i] = [(xoo - 
X3 3)(U)+2iIm( X03 )(Z A U)}/Tv[P +1 £(p)} and Trpp^] = 
[(Xn ~ X22)(U) - 2iRe(Xi2)(Z A U)]/Tr[P- 1 £(p)], where 
(Z A ), (U), and (Z A U) are all non-zero and already known. 
Therefore, we can obtain the four independent real parame- 
ters needed to calculate the coherence components X03 and 
X12, by simultaneously measuring, e.g., Z Z B and X A X B . 

In order to characterize the remaining coherence elements 
of x we make an appropriate change of basis in the prepa- 
ration of the two-qubit system. For characterizing xoi and 
X23, we can perform a Hadamard transformation on the first 
qubit, as H A \(j> c ) = a\+)\0) + /9|-)|1), where |±) = 
(|0) ± |l))/\/2. We then measure the stabilizer operator 
X A Z B , and a normalizer such as Z A X B . For characteriz- 
ing X02 and X3i, we prepare the system in the stabilizer state 
S A H A \(f> c ) = a\ +i}\0) +{3\ - and measure the stabi- 
lizer operator Y A Z B and a normalizer such as Z A X B , where 
S is the single-qubit phase gate, and |±i) = (|0) ± i \ 1))/ ^/2. 
Therefore, we can completely characterize the quantum dy- 
namical coherence with three BSM's overall. We note that the 
bases {|0), |1)}, {|±)} and {| ± i)} are mutually unbiased, 
i.e., the inner products of each pair of elements in these bases 
have the same magnitude J5| ■ 

A summary of the scheme for the case of a single qubit is 
presented in Fig. Q~]and Table. J] This table implies that the 
required resources in DCQD are as follows: (a) preparation of 
a maximally entangled state (for population characterization), 
(b) preparation of three other (nonmaximally) entangled states 
(for coherence characterization), and (c) a fixed Bell-state an- 
alyzer. We remark that a generalized DCQD algorithm for 



TABLE I: One possible set of input states and measurements for di- 
rect characterization of quantum dynamics (xij) f° r a single qubit, 
where |a| / \f3\ / 0, and {|0), {|±)}, {| ± t}} are eigenstates 
of the Pauli operators Z, X, and Y. 
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Measurement 


output 




Stabilizer 


Normalizer 




(|0)|0) + \1)\1))/V2 


Z A Z B ,X A X B 


N/A 


XOO, Xll, X22, X33 


a|0)|0)+/3|l)|l) 


Z A Z B 


X A X B 


X03, X12 


«I+>|0>+/3|-)|1) 


X A Z B 


Z A X B 


X01, X23 


a|+i)|0>+/3|-i>|l> 


Y A Z B 


Z A X B 


X02 , X13 
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TABLE II: Comparison of the required physical resources for char- 
acterizing a non-trace preserving CP quantum dynamical map on n 
qubits. N m and N m respectively denote the number of required input 
states and the number of noncommutative measurements for each in- 
put state, and N eX p = N- m N m is the number of required experimental 
configurations. 



Preparation 



Measurement 



Scheme 


dimCH) 


N ia 






SQPT 


2" 


4" 


4™ 


16" 


non-separable AAPT 




1 


4" + 1 


4' 1 + 1 


DCQD 


2"" 


-1" 


1 


4" 



qudits, with d being a power of a prime, is possible, and will 
be the subject of a future publication 0J. 

An additional feature of DCQD is that all the required en- 
semble measurements, for measuring the expectation values 
of the stabilizer and normalizer operators, can also be per- 
formed in a temporal sequence on the same pair of qubits with 
only one Bell-state generation. This is because at the end of 
each measurement, the output state is in fact in one of the four 
possible Bell states, which can be utilized as an input stabi- 
lizer state. 

For characterizing a quantum dynamical map on n qubits 
we need to perform a measurement corresponding to a ten- 
sor product of the required measurements for single qubits. 
An important example is a QIP unit with n qubits which has 
a 2™-dimensional Hilbert space (H). DCQD requires a total 
of 4™ experimental configurations for a complete characteri- 
zation of the dynamics. This is a quadratic advantage over 
SQPT and separable AAPT, which require a total of 16™ ex- 
perimental configurations. In general, the required state to- 
mography in AAPT could also be realized by non-separable 
(global) quantum measurements. These measurements can be 
performed either in the same Hilbert space, with 4" + 1 mea- 
surements, e.g., using a MUB measurement jj§], or in a larger 
Hilbert space, with a single generalized measurement j<§]. A 
comparison of the required physical resources is presented in 
Table HI1 A detailed resource cost analysis comparing DCQD 
to other QPT methods will be reported in a future publication 

i. 

Physical Realization. — For qubit systems, the resources re- 
quired in order to implement the DCQD algorithm are Bell 
state preparation and measurement, and single qubit rotations. 
These tasks play a central role in quantum information sci- 
ence, since they are prerequisites for quantum teleportation, 
quantum dense coding and quantum key distribution []]]. Be- 
cause of their importance, these tasks have been studied exten- 
sively and successfully implemented in a variety of different 
quantum systems, e.g., nuclear magnetic resonance (NMR) 
niOfl and trapped ions [11]. Thus, the DCQD algorithm is al- 
ready within experimental feasibility of essentially all systems 
which have been used to demonstrate QIP principles to date. 

Here, we propose a specific linear-optical implementation, 
realizable with present day technology - see Fig. [2] Using 
only linear-optical elements, at most 50% efficiency in dis- 
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FIG. 2: A schematic layout of a linear-optical realization of the 
DCQD algorithm. A pair of entangled photons, QHaHb) + 
e ltp |VaVb})/\/2, is generated via parametric down-conversion in 
a nonlinear crystal from a UV pump laser, where \H) and \ V) rep- 
resent horizontal and vertical polarization fl3l . The quarter and half 
waveplates are used to create non-maximally entangled states lfl3ll . 
and change the preparation and measurement bases. Two-photon in- 
terferometry occurs at the 50/50 beam splitters. This allows for dis- 
criminating two of the four Bell states deterministically, by utilizing 
polarizing beam splitters and photodetectors as analyzers. 



criminating among Bell states is possible [12]. The number of 
ensemble preparations in this specific setup is thus effectively 
increased by a factor of two over an implementation of the 
DCQD algorithm using an ideal Bell state analyzer. However, 
even in this optical realization, the number of experimental 
configurations is still reduced by a factor of 2™ (for character- 
izing a quantum dynamical process on n qubits) over SQPT 
and separable AAPT schemes. Alternatively, a small-scale 
and deterministic linear-optical implementation of DCQD can 
in principle be realized by a multirail representation of the 
qubits 11411 . 

To increase the efficiency of an all-optical Bell-state ana- 
lyzer, we have three different strategies: (i) introducing some 
nonlinearity, such as an optical switch 11511 . (ii) utilizing post- 
selected measurements, such as an optical CNOT jl<5], which 
has been demonstrated experimentally for filtering Bell states 
with a fidelity of about 79% lfl7ll . (iii) employing hyperentan- 
glement between pairs of photons for a complete and deter- 
ministic linear-optical Bell-state discrimination 

EEl. The 

latter method is based on the fact that the pairs of polarization- 
entangled photons generated by parametric down-conversion 
have intrinsic correlations in time-energy and momentum. 
These additional degrees of freedom can be exploited in order 
to distinguish between the subsets of Bell-states which other- 
wise cannot be discriminated by a standard Hong-Ou-Mandel 
interferometer. Note that employing the time-energy correla- 
tions in the context of implementing the DCQD algorithm is 
feasible if the quantum dynamical map acts only on the polar- 
ization degrees of freedom. 

Partial Characterization of Dynamics. — DCQD can be ap- 
plied, by construction, to the task of partial characterization 
of quantum dynamics, where we cannot afford or do not need 
a full characterization of the system, or when we have some 
a priori knowledge about the dynamics. In particular, we 
can substantially reduce the number of measurements, when 
we are interested in estimating the coherence elements of the 
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superoperator for only specific subsets of the operator ba- 
sis and/or subsystems of interest. E.g., we need to perform 
a single ensemble measurement if we are required to iden- 
tify only the coherence elements X03 an d X12 of a particular 
qubit. Here, we present an example of such a task. Specif- 
ically, we demonstrate that the DCQD algorithm enables the 
simultaneous determination of coarse-grained physical quan- 
tities, such as the longitudinal relaxation time T\ and the trans- 
verse relaxation (or dephasing) time T2. Assume that we pre- 
pare a two-qubit system in the non-maximally entangled state 
\4>c) = a 1 00) + /3|11). Then we subject qubit A to an am- 
plitude damping process for duration t\, followed by a phase 
damping process, for duration t2- The elements of the final 
density matrix then read (0|p/|0) = 1 - exp(— — H 2 ) 

and (0|/0/|l) = exp(— )a*(3, where ^7 = || + In or- 
der to determine T\, we measure the eigenvalues of the stabi- 
lizer operator Z A Z B . We obtain either +1 or —1, correspond- 
ing to the projective measurement P+i, or P-i. The proba- 
bilities of either of these outcomes, e.g., Tr(P-ipf) are re- 
lated to Ti through the relation i = -± ln(l - gg|g^g4 ). 
In order to obtain information about T%, we measure the ex- 
pectation value of any normalizer of the input state, such as 

X A X B , yielding jf? = -2 In ^*x*x>?'p) • Since the °P" 
erators Z A Z B and X A X B commute, we can measure them 
simultaneously. Therefore we can find both T\ and T 2 in a 
Bell-state measurement. 

Outlook. — One can combine the DCQD algorithm with the 
method of maximum-likelihood estimation 12011 . in order to 
minimize the statistical errors in each experimental configu- 
ration. Moreover, a new scheme for continuous characteri- 
zation of quantum dynamics can be introduced, by utilizing 
weak measurements for the required quantum error detections 
in DCQD J2H. 

For quantum systems with controllable two-body interac- 
tions (e.g., trapped-ion and NMR systems), DCQD could have 
near-term experimental applications for complete verification 
of small QIP units. For example, DCQD, reduces the number 
of required experimental configurations for systems of 3 or 4 
physical qubits from 5 x 10 3 and 6.5 x 10 4 (in SQPT) to 64 
and 256, respectively. A similar scale-up can only be achieved 
by utilizing non-separable AAPT methods. Complete char- 
acterization of such dynamics would be essential for verifi- 
cation of quantum key distribution procedures, teleportation 
units, quantum repeaters, and more generally, in any situation 
in quantum physics where a few qubits have a common local 
bath and interact with each other. Furthermore, as demon- 
strated here, DCQD is inherently suited to extract partial in- 
formation about quantum dynamics. Several other examples 
of such applications have been demonstrated. Specifically, it 



has been shown that DCQD can be used for realization of gen- 
eralized quantum dense coding |8y. Moreover, DCQD can 
be efficiently applied to (single- and two-qubit) Hamiltonian 
identification tasks |@]. Finally, the general techniques devel- 
oped here could be further utilized for closed-loop control of 
open quantum systems. 
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